Generalized Parton Distributions for the Proton in Position Space 

Non-Zero Skewness 



R. Manohar", A. Mukherjee", D. Chakrabarti^ 

Department of Physics, 
Indian Institute of Technology Bombay, 
Powai, Mumbai 400076, India 
^ Department of Physics, 
Indian Institute of Technology Kanpur, 
Kanpur 208016, India 
(Dated: December 14, 2010) 

Abstract 

We investigate the generalized parton distributions (GPDs) for u and d quarks in a proton in 
transverse and longitudinal position space using a recent phenomenological parametrization. We take 
nonzero skewness C, and consider the region x > Q. Impact parameter space representation of the GPD 
E is found to depend sharply on the parameters used within the model, in particular in the low x 
region. In longitudinal position space a diffraction pattern is observed, as seen before in several other 
models. 
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I. INTRODUCTION 



The generalized paton distributions (GPDs) have gained a lot of theoretical and experimen- 
tal interest recently. Unlike the ordinary parton distributions (pdfs) which at a given scale 
depend only on the longitudinal momentum fraction x of the parton, GPDs are functions of 
three variables, x, ( and t where the so-called skewness ( gives the longitudinal momentum 
transfer and —t is the square of the momentum transfer in the process. The GPDs give inter- 
esting information about the spin and orbital angular momentum of the quarks and gluons in 
the nucleon. They are experimentally accessed through the overlap of deeply virtual Compton 
scattering (DVCS) and Bethe-Heitler (BH) process as well as exclusive vector meson produc- 
tion [1]. The real and imaginary parts of the Compton amplitude give information on GPDs 
in different kinematical regions. Several experiments worldwide, for example, at DESY HERA 
collider, by the HI [2, 3| and ZEUS 0, 0] collaboration and HERMES jsj fixed target experi- 
ments have finished taking data on DVCS. Experiments are also being done at JLAB Hall A 
and B COMPASS at CERN has programs to access GPDs through muon beams [8]. A 

review of the different models and their status with respect to the experimental results can be 
ft 

found in [9|. As the GPDs involve a momentum transfer (off- forwardness) , they do not have 
probabilistic interpretation, unlike ordinary parton distributions (pdfs). However, it has been 
shown that when the momentum transfer is purely in the transverse direction, if one performs 
a Fourier Transform (FT) with respect to the transverse momentum transfer A^, one gets the 
so-called impact parameter dependent parton distributions (ipdpdfs), originally proposed by 



Soper in the context of nucleon form factor 
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Ipdpdfs tell us how the quarks of a given 



longitudinal momentum are distributed in the transverse position or impact parameter space. 
These obey certain positivity conditions and unlike the GPDs, have probabilistic interpreta- 
tion. As transverse boosts are non-relativistic Galilean boosts in light-front formalism, there 
is no relativistic correction to this interpretation, even though we are considering relativistic 
systems. Ipdpdfs are defined in a proton state with a sharp plus momentum p'^ and localized 
in the transverse plane such that the transverse center of momentum R± = (normally, one 
should work with a wave packet state which is very localized in transverse position space, in 
order to avoid the state to be normalized to a delta function 1^). These give simultaneous 
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information about the longitudinal momentum fraction x and the transverse distance h of the 
parton from the center of the proton and thus give a new insight to the internal structure of 
the proton. A wave packet state which is transversely polarized is shifted sideways in impact 
parameter space. This shift is determined by the GPD E. An interesting interpretation of Ji's 
sum rule is given in |pL^] in terms of ipdpdfs, with E related to the orbital angular momentum 
carried by the quarks. Interesting connections between the ipdpdfs and transverse momentum 
dependent pdfs have been obtained in various models. However such relations cannot exist in 
a model independent way jl^ . 

Since GPDs depend on a sharp x, the Heisenberg uncertainty relation restricts the longitu- 
dinal position space interpretation of GPDs themselves. It has, however, been shown in 15 1 
that one can define a quantum mechanical Wigner distribution for the relativistic quarks and 
gluons inside the proton. Integrating over k~ and fc"*", one obtains a four dimensional quantum 
distribution which is a function of r and where r is the quark phase space position defined 
in the rest frame of the proton. These distributions are related to the FT of GPDs in the 
same frame. This gives a 3D position space picture of the GPDs and of the proton, within 



the limitations mentioned above. In 



161], a parametrizations of generalized parton correlation 



functions have been done for a spin 1/2 target. These reduce to GPDs when integrated over 
k± and to TMDs when the momentum transfer is zero. 

In another work, the real and imaginary parts of the DVCS amplitude are expressed in 
longitudinal position space by introducing a longitudinal impact parameter a conjugate to 
the skewness C,. Taking a field theory inspired simple relativistic spin 1/2 system, namely 
for an electron dressed with a photon in QED, it was found that the DVCS amplitude show a 
diffraction-like pattern in longitudinal position space [l3] • Since Lorentz boosts are kinematical 
in the front form, the correlation defined in the 3 D position space 6_l and a is frame independent. 
Similar diffraction pattern was observed in a holographic model for the meson. In this work, we 
investigate a recent parametrization [l^ of the GPDs in position space for non-zero C,. GPDs 
for zero skewness have been parametrized in a similar way in 19|]. At the input scale they are 
parametrized by a spectator model term multiplied by a Regge motivated term. The parameters 
were fitted by fitting the forward pdfs and form factors. For non-zero the GPDs have to satisfy 
an additional constraint, namely polynomiality. In certain models, for example using the overlap 
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of light- front wave functions (LFWFs), it is very difficult to obtain a suitable parametrization 
of the higher Fock components of the wave function in order to get the polynomiality of GPDs. 
Polynomiality is satisfied by construction only if one considers the LFWFs of simple spin 1/2 
objects like a dressed quark or a dressed electron in perturbation theory instead of the proton 



2l|. A recent fit to the DVCS data at small Bjorken x from HI and ZEUS was done in 



22], using the conformal Mellin-Barnes representation of the DVCS amplitude. However, to 
get the GPDs one has to do an inverse Mellin transform, and a knowledge of all moments are 
required for that. In 18j a functional form of GPDs in the DGLAP region x > ( has been 
obtained by generalizing a parametrization for zero skewness. In the ERBL region x < (, 
a quark and an antiquark pair emerges from the nucleon and undergoes the electromagnetic 
interaction. Lattice moments were used as constraints and weighted average of GPDs around 
definite values of x were constructed using Bernstein polynomials. However, so far only a few 
of the lattice moments are known and this method has large theoretical uncertainties. 

I. a p.ev,ous pubHcation Q we used the GPD pa..a.e.izati„n of Q at .e.o s.ew„e.. to 
calculate the distributions of partons in the transverse impact parameter space. In this work, 
we study the GPDs in transverse and longitudinal position space using the phenomenological 
parametrization in \18l for nonzero skewness. The plan of the paper is as follows. After 
presenting the model used in section II, we discuss numerical results for the GPDs in transverse 
and longitudinal position space respectively in section HI. Summary and conclusions are given 
in section IV. 



II. PARAMETRIZATION OF THE GPDS 



We consider the parametrization in 




for the GPDs for nonzero skewness (C 7^ 0) : 



Set I 

H\x,C,t) = Gj;J,(x,C,t)a;-"'~''i(^-")'' (1) 
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Set II 



E'\x, C, t) = C, (4) 

X 

All parameters except for pi and p2 are flavor dependent. The function G has the same form 
for both parametrizations, I and II: 



X 



GUx.Ct)= AT— (5) 



where 



D{x,k^) = e-m^, (6) 



and 



A;2 = - -^M^ - (7) 

1 - X 1 - X ^ ^ 

,'^ = ,M^__4_Mj-(k,-i^A)=i^, (8) 



When the skewness is nonzero, the total momentum transfer square is modified to 



^ = = + (1 - e)D' = -t^^n + (1 - e)D' (10) 



where D = PV(1 - ^ - ^/(l + and ^ = C/{2 - ()■ D reduces to A at C = 0. 

Here x is the fraction of the light cone momentum carried by the active quark, k being its 
momentum. The mass parameters are m, the struck quark mass, and M, the proton mass. 
The normalization factor includes the nuclcon- quark- diquark coupling, and it is set to A/" = 1 
GeV^. Here we consider the DGLAP region < x < 1. The dominating contribution in this 
kinematical region comes from the process where a quark from the proton with momentum 
fraction x is struck by the incident photon and again reabsorbed by the proton. The above 



phenomenologically motivated parametrization of the GPDs H{x,(,t) and E{x,(,t) at zero 
skewness ( was done in [l9| using a spectator model calculation at the low input scale. The 
spectator model has been used for its simplicity and for the fact that it is flexible enough 
to predict the main features of a number of distribution and fragmentation functions in the 
intermediate and large x region. The spectator mass is chosen to be different for different 
quark flavor GPDs. However, similar to the case of pdfs, the spectator model is not able to 
reproduce quantitatively the small x behaviour of the GPDs. So a 'Regge-type' term has been 
considered multiplying the spectator model function G^f^. Extension to nonzero ( was done 
in [l8|. The parameters are listed in IiqI for both the sets. The parameters M^, A"^ and a^, 
q = u,d, obtained at an initial scale Qq {Qq = 0.094 GeV^), and they are the same for both 
Sets I and II, in Set I they are by definition the same for the functions H and E (see Eqs. f[TP|) ). 
The parameters (32, pi and p2, in Set I, and all parameters defining E in Set II (Eq. (jl]), 
were fitted to the nucleon electric and magnetic form factors, with the values of M|, A^, and 
fixed. The input scale Qq = 0.094GeV^ is obtained as a parameter in the model. The low 
value of Qq results from the requirement that only valence quarks contribute in the momentum 
sum rule. The GPD H is constrained in the forward limit by the pdf data. As the GPD E is 
unconstrained by the data on forward pdfs, two different sets of parameters were used in the 
fit, these are denoted by set I and II. In set II an additional normalization condition was used 

dxEg{x,t = 0) = k'^ (11) 



where and k"^ are the u and d quark contributions to the nucleon anomalous magnetic 
moment. Although and H'^ have similar behaviour in both the sets, E"^ and E'^ behave 
differently even for C = 0. In the forward limit for H, Alekhin leading order pdf sets were 
used. Additional constraints for nonzero ( were obtained from lattice moments. 



III. PARTON DISTRIBUTIONS IN IMPACT PARAMETER SPACE 

Parton distribution in impact parameter space is defined as : 

q{xX,b) = ^1 d'A^e-'^--'-H{xX,t) 
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e(x, C,b) = ^l d'A^e-'^^-'-Eix, C, t). (12) 

Here 6 =| 6^ | is the transverse impact parameter which is a measure of the transverse 
distance between the struck parton and the center of momentum of the hadron. b is defined 
such a way that ^i^i = where the sum is over the number of partons. The relative distance 
between the struck parton and the spectator system provides an estimate of the size of the 
system as a whole. The above picture was proposed in the limit of zero skewness ( in In 
most experiments ( is nonzero, and it is of interest to investigate the GPDs in b± space with 
nonzero (. 

As described in the introduction the ipdpdfs describe the probability of finding a parton 
of definite momentum fraction a; at a distance b± from the center of the proton. When the 
skewness ( is nonzero, the transverse location of the proton itself is different before and after 
the scattering. This transverse shift does not depend on x but on the skewness ( and b = \ b± \. 
The information on the transverse shift is not washed out even if the GPDs are integrated over 
X in the DVCS amplitude. In the DGLAP region x > (, the impact parameter b± gives the 
location where the quark is pulled out and put back to the nucleon. In the ERBL region ( < x, 
b± denotes the transverse location of the quark-antiquark pair inside the nucleon. For a single 
fermion, the impact parameter dependent pdf would be a delta function. The smearing in bj_ 
space is due to the multiparticle correlation. In Fig. 1 we plot eu{x,(,b) and ed{x,(,b) as a 
function of b for a fixed value of x and different values of (. Substantial difference is seen in the 
impact parameter space for the two sets of parametrization, set I and set II. In set I both and 
Cd become more sharply peaked at 6 = when ( is smaller. It is to be noted that the change 
in e{x, b) with b is related to the deformation of the parton distribution in impact parameter 
space for transversely polarized nucleon. For non-zero this probes the deformation when the 
nucleon has a transverse shift before and after scattering. In set II for e„, the peak increases 
as C increases, in contrast to set I, whereas for e^, the peak decreases for increasing but does 
not become broader as in set I. 

In Fig. 2, we have plotted C, b) as a function of b for fixed x = 0.6 and different values of 
(. As the two parametrizations set I and II are not much different for the GPD H{x,C,,t) not 
much difference is seen in the impact parameter space. In all cases here the peak in C, b) for 
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fixed X decreases with increase of (. For d quark distributions, tlie peak also becomes broader 
as ( increases, wliicli is tlie same befiaviour as seen in e(a;,C,6). As ( increases for fixed x, 
the smearing in b space becomes broader, wliicli means tliat as tlie momentum transfer in the 
longitudinal direction increases, the active quark is more likely to be pulled out at a larger b. 

In Figs. 3 and 4 we plotted the impact parameter dependent distributions e{x,(, b) and 
q{x, (, b) for fixed ( and b values as functions of x. Substantial difference is observed in e{x, (, b) 
between set I and II, in particular in the low x region. For the peak shifts towards smaller 
X values {x > ( in our calculation). However for — the peak shifts to higher x values in set 
II. The qualitative behaviour of g„ and qa are the same as functions of x. ^ ^ 

The boost invariant longitudinal impact parameter a was first introduced in il7|| and it was 
shown that DVCS amplitude shows interesting diffraction pattern in longitudinal impact pa- 
rameter space. GPDs also when expressed in term of a exhibit the similar diffraction pattern 



251]. The boost invariant longitudinal impact parameter conjugate to the longitudinal momen- 



tum transfer is defined as a = |6 P^. So, the GPDs in longitudinal position space is given 
by: 



q{x,a,t) = ^ ! ' dCe<''^'-'^H{xX.t) 
= ^1^' dCe'^''H{xX,t) 



) 

e(x,a,t) = ^ dCe^'-^'-l'EixX.t) 

= ^ r dCe'^''E{xX,t) (13) 

^TT Jo 



Since we are concentrating only in the region ^ < x < 1, the upper limit of C, integration Cj is 
given by Cmax if x is larger then Cmax, otherwise by x if x is smaller than C,max where Cmax is the 
maximum value of C, allowed for a fixed —t: 



i-t) ( \ 4M2 . ^ ^ 



2Mny (-t) 

In Figs. 5 and 6 we have plotted the GPDs in longitudinal position space a . We restrict 
ourselves to the DGLAP region. The GPDs show diffraction pattern in a space, similar to that 

n 

observed for a dressed electron in QED or in a holographic model for the meson [17J. There is 



a primary maximum followed by a series of secondary maxima. The positions of the minima 
are the same for both u and d quark GPDs and they do not depend whether we are plotting 
q{x,a,b) or e{x,a,b). The feature was also observed for a dressed electron state in [l7|. The 
positions of the minima are characteristics of the finite Fourier transform and independent of 
the GPD used. As —t increases, the positions of the first minima move in to smaller values of 
(J. As seen in Fig. 5, for set II, the peaks in eu{x,a,t) are more sharp compared to set I. On 
the other hand, for ed{x, a, t) the magnitude of the peak is more affected as —t increases in set 
II compared to set I. However for qu{x, a, t) and qd{x, a, t) there is not much difference between 
the two parametrization sets, as expected. In all cases, the magnitude of the peak decreases as 
—t increases and the first minima move in. This effect has already been observed for a dressed 
electron state and in a holographic model for the meson as well as for chiral odd GPDs 



25| . In 171] a relation between the position of the first minima and the momentum transfer 
squared —t has been derived in analogy with diffraction in optics. The fact that such pattern 
with the same qualitative behaviour is observed here in a phenomenological parametrization 
of GPDs shows that it is not an artefact of the cuts and constraints used in the field theory 



based model. In j26| pion distribution amplitudes (DAs) have been investigated in longitudinal 
position space (loffe time) by expressing the DAs in an expansion over Gegenbauer polynomials. 
Oscillation of the DA has been observed in longitudinal position space showing two humps in 
certain models. It was concluded that this is an artefact of truncating the conformal spin partial 
wave expansion. In order to further investigate the diffraction pattern that we observe for the 
GPDs, in Fig. [7] we plot the GPD H in longitudinal position space using a model similar to 
that in . For the dotted and long dashed curves we chose a factorized ansatz of the form 

H{xX,t) = Nu{x)Fi{t) (15) 

where is the normalization constant, we took u{x) = x{l~x) and for the form factor Fi{t) we 
took the standard dipole form. The FT shows diffraction pattern with the general features the 
same as observed in other models. In the solid and short-dashed curves, we used a somewhat 
different ansatz, namely 

H{xX,t)=Afu{x,OFi{t) (16) 
where we took the same dipole form factor, and u{xX) = {x — C)^(l ^ ^)^; being the 
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normalization constant. The normalization constants have been chosen to fit the two sets of 
curves in the same plot. There is no diffraction pattern in the second model; which clearly 
shows that the pattern is not only due to the finite range of the ( integration but depends also 
on the X, ( and t interplay in the GPD parametrization used. So an experimental observation of 
such diffraction pattern can help to constrain GPD parametrizations. However in order to get 
the full Lorentz invariant picture in longitudinal position space one has to consider the other 



kinematical region x < ( as well; as explained in 



ITj . Also the positions of the minima and 



their shift with increase of —t are model independent features of the Fourier transform. 



IV. CONCLUSION 

In this work we investigated the GPD for u and d quark distributions in the proton using a 
recent parametrization [18]. For nonzero ( we worked in the DGLAP region x > (. Taking a 
Fourier transform with respect to the transverse momentum transfer A_l we obtained the parton 
distributions in impact parameter space. For nonzero ( these probe the parton distributions 
when the initial proton is shifted from the final proton in the transverse plane. When the 
proton is transversely polarized, the parton distributions in the tranverse plane is distorted, 
this distortion is related to the GPD E{x,(,t) and to the orbital angular momentum of the 
quarks. We showed the x and bj_ dependence of the ipdpdfs for various ( values. We introduced 
a boost invariant longitudinal impact parameter a conjugate to C. Both the GPDs H and E in 
cr space show diffraction pattern as seen before in some other models. We did some comparative 
study using different toy models and showed that although the general features of this pattern 
are independent of specific models but the appearance of the diffraction pattern depends not 
only on the finiteness of the ( integration but also on the interplay of the x, ( and t dependence 
of the GPDs. 
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FIG. 1: (Color online) Plots of (a) eu{x, C,b) vs b =\ b± \ for fixed values of x and C and for parameters 
as in set I (b) same as in (a) but for d quark, (c) same as in (a) but the parameters are as in set II, 
(d) same as in (b) but the parameters are as in set II. b is in GeV~^. 
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(c) b (d) b 

FIG. 2: (Color online) Plots of (a) ^^(a;, C, vs 6 = | 6_l | for fixed values of x and C, and for parameters 
as in set I, (b) same as in (a) but for d quark, (c) same as in (a) but the parameters are as in set II, 
(d) same as in (b) but for parameters as in set II. h is in GeV~^. 
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FIG. 3: (Color online) Plots of (a) eu{x, C, b) vs x for fixed values of b and C and for parameters as in 
set I, (b) same as in (a) but for d quark, (c) same as in (a) but for parameters as in set II, (d) same 
as in (b) but for parameters as in set II. b is in GeV~^. 
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FIG. 4: (Color online) Plots of (a) Quix, C, b) vs x for fixed values of b and ( and for parameters as in 
set I, (b) same as in (a) but for d quark, (c) same as in (a) but for parameters as in set II, (d) same 
as in (b) but for parameters as in set II. b is in GeV""*^. 
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